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Abstract 

The aim of the paper is to propose a bounded-error quantum polyno¬ 
mial time (BQP) algorithm for the max-bisection and the min-bisection 
problems. The max-bisection and the min-bisection problems are funda¬ 
mental NP-hard problems. Given a graph with even number of vertices, 
the aim of the max-bisection problem is to divide the vertices into two 
subsets of the same size to maximize the number of edges between the 
two subsets, while the aim of the min-bisection problem is to minimize 
the number of edges between the two subsets. The proposed algorithm 
runs in 0{rn?) for a graph with m edges and in the worst case runs in 
O(n^) for a dense graph with n vertices. The proposed algorithm targets a 
general graph by representing both problems as Boolean constraint satis¬ 
faction problems where the set of satisfied constraints are simultaneously 
maximized/minimized using a novel iterative partial negation and partial 
measurement technique. The algorithm is shown to achieve an arbitrary 
high probability of success of 1 — e for small e > 0 using a polynomial 
space resources. 

Keywords: Quantum Algorithm, Graph bisection. Max-bisection, Min- 
bisection, Amplitude Amplification, BQP, NP-hard. 


1 Introduction 

Given an undirected graph G = (R E) with a set V of even number of vertices 
and a set E of unweighted edges. Two graph bisection problems will be consid¬ 
ered in the paper, the max-bisection problem and the min-bisection problem. 
The goal of the max-bisection problem is to divide V into two subsets A and 
B of the same size so as to maximize the number of edges between A and B, 
while the goal of the min-bisection problem is to minimize the number of edges 
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between A and B. In theory, both bisection problems are NP-hard for general 
graphs [Til[5]. 

These classical combinatorial optimization problems are special cases of 
graph partitioning m- The graph partitioning has many applications, for ex¬ 
ample, divide-and-conquer algorithms [24], compiler optimization |21j . VLSI cir¬ 
cuit layout |5], load balancing image processing [30] . computer vision [22] . 
distributed computing [25] . and route planning [Tj. In practice, there are many 
general-purpose heuristics for graph partitioning, e.g. (TS] |3I1 131] that han¬ 
dle particular graph classes, such as UniElEl]- There are also many practical 
exact algorithms for graph bisection that use the branch-and-bound approach 
[23l |8]. These approaches make expensive usage of time and space to obtain 
lower bounds [UHKHiis]. 

On conventional computers, approximation algorithms have gained much 
attention to tackle the max-bisection and the min-bisection problems. The 
max-bisection problem, for example, has an approximation ratio of 0.7028 due 
to |10j which is known to be the best approximation ratio for a long time by in¬ 
troducing the RPR^ rounding technique into semidefinite programming (SDP) 
relaxation. In |12] . a poly-time algorithm is proposed that, given a graph ad¬ 
mitting a bisection cutting a fraction 1 — e of edges, finds a bisection cutting 
an (1 — g{e)) fraction of edges where g{£) —>■ 0 as e —>■ 0. A 0.85-approximation 
algorithm for the max-bisection is obtained in [28]. In [33], the SDP relaxation 
and the RPR^ technique of m have been used to obtain a performance curve 
as a function of the ratio of the optimal SDP value over the total weight through 
finer analysis under the assumption of convexity of the RPR^ function. For the 
min-bisection problem, the best known approximation ratio is 0{log n) [27] with 
some limited graph classes have known polynomial-time solutions such as grids 
without holes m and graphs with bounded tree width [20] . 

The aim of the paper is to propose an algorithm that represents the two 
bisection problems as Boolean constraint satisfaction problems where the set of 
edges are represented as set of constraints. The algorithm prepares a superposi¬ 
tion of all possible graph bisections using an amplitude amplification technique 
then evaluates the set of constraints for all possible bisections simultaneously 
and then amplifies the amplitudes of the best bisections that achieve the max¬ 
imum/minimum satisfaction to the set of constraints using a novel amplitude 
amplification technique that applies an iterative partial negation and partial 
measurement. The proposed algorithm targets a general graph where it runs in 
0{nn?) for a graph with m edges and in the worst case runs in O(n^) for a dense 
graph with number edges close to m = ^ vertices to achieve an 

arbitrary high probability of success of 1 — e for small e > 0 using a polynomial 
space resources. 

The paper is organized as follows; Section 2 shows the data structure used to 
represent a graph bisection problem as a Boolean constraint satisfaction prob¬ 
lem. Section 3 presents the proposed algorithm with analysis on time and space 
requirements. Section 4 concludes the paper. 


2 Data Structures and Graph Representation 

Several optimization problems, such as the max-bisection and the min-bisection 
problems, can be formulated as Boolean constraint satisfaction problems 013] 
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where a feasible solution is a solution with as many variables set to 0 as variables 
set to 1, i.e. balanced assignment, as follows: for a graph G with n vertices 
and m edges, consider n Boolean variables vq,. .. ,Vn-i and m constraints by 
associating with each edge (a, b) € E the constraint c; = Va (B Vb, with I = 
0,1,..., m — 1, then the max-bisection is the problem that consists of finding 
a balanced assignment to maximize the number of constraints equal to logic-1 
from the m constraints, while the min-bisection is the problem that consists of 
finding a balanced assignment to maximize the number of constraints equal to 
logic-0 from the m constraints, such that if a Boolean variable is set to 0 then 
the associated vertex belongs to the first partition and if a Boolean variable is 
set to 1 then the associated vertex belongs to the second partition. 



Figure 1: (a) A random graph with 8 vertices and 12 edges, (b) A max-bisection 
instance for the graph in (a) with 10 edges connecting the two subsets, and (c) 
(b) A min-bisection instance for the graph in (a) with 3 edges connecting the 
two subsets. 


For example, consider the graph G shown in Figure [IJa). Let G = {V,E), 
where, 

ld= {0,1, 2, 3,4, 5, 6, 7}, 

A = 1(0,1), (0,2), (0,3), 

(1,2), (1,7), (2, 3), (1) 

(3,4), (3,6), (4,5), 

(4,6), (5,7), (6,7)}. 

Assume that each vertex a S 17 is associated with a Boolean variable Va, 
then the set of vertices V can be represented as a vector X of Boolean variables 
as follows, 

X = {vo,vi,V2,V3,V4,V5,ve,vr), (2) 
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and if each edge (a, b) € E is associated with a constraint ci = Va(B Vb then the 
set of edges E can be represented as a vector Z of constraints as follows, 


z = (co, Cl, C2, C3, C4, C5, C6, C7, Cg, Cg, Cig, Cn), 


(3) 


such that, 

Co = {vo © Vi),Ci = {vo © V 2 ),C 2 = {vq © V 3 ), 

C3 = (vi © V 2 ), C4 = {vi © Vr), C5 = {V 2 © V 3 ), 

cq = (C3 © Vi), C7 = {V 3 © Ce), Cg = (C4 © r;5), 

eg = {V 4 © Ve), Cio = (C5 © 'ey), Cii = (vq © V 7 ) 


In general, a bisection Gp for the graph G can be represented as Gp = 
(x, z(x)) such that each vector x G {0,1}" of variable assignments is associated 
with a vector z(x) € {0,1}"* of constraints evaluated as functions of the variable 
assignment x. In the max-bisection and the min-bisection problems, the vector 

/ n 

X of variable assignments are restricted to be balanced so there are M = „ 

\ 2 

possible variable assignments among the N = 2^ possible variable assignments, 
and the solution of the max-bisection problem is to find the variable assign¬ 
ment that is associated with a vector of constraints that contains the maximum 
number of I’s, and the solution for the min-bisection problem is to find the 
variable assignment that is associated with a vector of constraints that contains 
the maximum number of O’s. For example, for the graph G shown in Figure 
[TJa), a max-bisection for G is ((0,1, 0,1, 0,1,1,0), (1,0,1,1,1,1,1, 0,1,1,1,1)) 
with 10 edges connecting the two partitions as shown in Figure db), and a 
min-bisection for G is ((0, 0, 0,0,1,1,1,1), (0,0, 0,0,1,0,1,1, 0,0,0,0)) with 3 
edges connecting the two partitions as shown in Figure [TJc). It is important 
to notice that a variable assignment x = (0,1, 0,1, 0,1,1,0) is equivalent to 
X = (1, 0,1, 0,1, 0,0,1), where x is the bit-wise negation of x. 


3 The Algorithm 

Given a Graph G with n vertices and m edges. The proposed algorithm is di¬ 
vided into three stages, the first stage prepares a superposition of all balanced 
assignments for the n variables. The second stage evaluates the m constraints 
associated with the m edges for every balanced assignment and stores the values 
of constraints in constraint vectors entangled with the corresponding balanced 
assignments in the superposition. The third stage amplifies the constraint vec¬ 
tor with maximum (minimum) number of satisfied constraints using a partial 
negation and iterative measurement technique. 


3.1 Balanced Assignments Preparation 

To prepare a superposition of all balanced assignments of n qubits, the proposed 
algorithm can use any amplitude amplification technique, e.g. [13 [Ml isg. 
An extra step should be added after the amplitude amplification to create an 
entanglement between the matched items and an auxiliary qubit |axi), so that 
the correctness of the items in the superposition can be verified by applying 
measurement on |axi) without having to examine the superposition itself. So, 
if \axi) = 11) at the end of this stage, then the superposition contains the correct 
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Figure 2: A quantum circuit for the proposed algorithm. 


items, i.e. the balanced assignments, otherwise, repeat the preparation stage 
until |aa;i) = |1). This is useful for not having to proceed to the next stages 
until the preparation stage succeeds. 

The preparation stage to have a superposition of all balanced assignments 
of n qubits will use the amplitude amplification technique shown in |35) since it 
achieves the highest known probability of success using fixed operators and it 
can be summarized as follows, prepare a superposition of 2" states by initializing 
n qubits to state |0) and apply iJ®” on the n qubits 


|5-o) = (H^^) |0)®" 

iV-1 


(5) 


where H is the Hadamard gate, and N = 2”. Assume that the system ITq) is 
re-written as follows, 


N-l N-l 

^ ^ ^ (6) 

j—0, 3—0, ^ ^ 

JGXt j^Xp 

where Xt is the set of all balanced assignments of n bits and Xp is the set 
of all unbalanced assignments. Let M = ^ ^ ^ be the number of balanced 

assignments among the 2" possible assignments, sin(0) = y/M/N and 0 < 0 < 
7r/2, then the system can be re-written as follows, 

l^'o) = sin(6») I'f/;!)-kcos(6») IV’o), (7) 

where \tjji) = |r) represents the balanced assignments subspace and \ijjo) repre¬ 
sents the unbalanced assignments subspace. 

Let D = kFi?o ((/-) W^Rr {<!)), i?o (</>) =/-(1 - |0) (0|, Rr (</>)=/-(!- 

e®'^) |t) (r|, where W = is the Walsh-Hadamard transform [T^. Iterate the 
operator D on jTo) for q times to get, 

|vI/i)=i?«|vI/o) = a,|V'i)+MV'o), (8) 
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such that, 


(9) 


a, = sin(0) {y) + {y)) , 


b, = cos(0)e*(«-i)'^ (U, {y) + C/,_i (y)), (10) 

where y = cos((5), cos (J) = 2 sin^(0) sin^(^) — 1, 0 < 0 < 7r/2, and Uq is the 
Chebyshev polynomial of the second kind [55] defined as follows, 


Uq {y) = 


sin((g + 1) (5) 


Setting (j) = 6.02193 


sin (S) 

1.91687r, M = 


( 11 ) 


, iV = 2" and, q = 


sm{9) 


J. 


then \aq\ > 0.9975 


] . The upper bound for the required number of iterations 


q to reach the maximum probability of success is. 


9 = 


sin(0) J 


< 1.91687r 



( 12 ) 


and using Stirling’s approximation. 



(13) 


then, the upper bound for required number of iterations q to prepare the super¬ 
position of all balanced assignments is. 


q Ri l.9168y —n = O (v^) • (14) 

It is required to preserve the states in lipi) for further processing in the next 
stage. This can be done by adding an auxiliary qubit \axi) initialized to state 
|0) and have the states of the balanced assignments entangled with |axi) = |1), 
so that, the correctness of the items in the superposition can be verified by 
applying measurement on |aa:i) without having to examine the superposition 
itself. So, if \axi) = |1), then the superposition contains the balanced assign¬ 
ments, otherwise, repeat the preparation stage until \ax\) = |1). This is useful 
to be able to proceed to the next stage when the preparation stage succeeds. 
To prepare the entanglement, let 

| 4 ' 2 ) = | 0 ) , , 

= aq IV'i) 0 |0) -b bq I'ipo) ® |0) , ^ 

and apply a quantum Boolean operator Uf on |4'2), where Uf is defined as 
follows. 


Uf |a;,0) 


|a;,0) ,if |x) e [V'o), 
lx, 1) ,if |x) G IV'i), 


(16) 


and / : {0,1}" —>■ {0,1} is an n inputs single output Boolean function that 
evaluates to True for any x G Xt and evaluates to False for any x £ Xp, then. 


\'^s)=Uf\^2) 

= aq IV'i) 0\1) + bq IV'o) ® |0). 


(17) 
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Apply measurement Mi on the auxiliary qubit \axi) as shown in Figure [21 


The probability of finding \axi) = |1) is. 


Pr{Mi = 1) = |a,|^ > 0.9975, 

(18) 

and the system will collapse to. 



(19) 


3.2 Evaluation of Constraints 


There are M states in the superposition 




(Mi = l) 


each state has an amplitude 


then let |\E' 4 ) be the system after the balanced assignment preparation stage 
as follows, 


M-l 

|vl/4)=a5^ Ixfc), (20) 

k=0 

where |aa;i) is dropped from the system for simplicity and a = For a 

graph G with n vertices and m edges, every edge (a, b) connecting vertcies 
a,b € y is associated with a constraint c/ = Va (B Vb, where Va and Vb are 
the corresponding qubits for vertices a and b in |'I' 4 ) respectively such that 
0<Z<m, 0<m< , 0<a,6<n — 1 and a ^ b, where is the 

maximum number of edges in a graph with n vertices. 

To evaluate the m constraints associated with the edges, add m qubits ini¬ 
tialized to state |0), 


|^' 5 ) = |^' 4 )®| 0 )®’” 
k=0 

For every constraint c; = Va(Bvb, apply two Cont_ax gates, Cont_ax{va,ci) 
and Cont_ax{vb,ci), so that |c/) = |ua0Uh). The collection of all Cont^ax 
gates applied to evaluate the m constraints is denoted Cy in Figure |2l then the 
system is transformed to, 


M-l 

1^6) = a ^ (Ixfc)® |cgcj...ct_i)), (22) 

fe=0 

where ax is the Pauli-X gate which is the quantum equivalent to the NOT gate. 
It can be seen as a rotation of the Bloch Sphere around the X-axis by tt radians 
as follows. 


<yx 


0 1 
1 0 ’ 


(23) 


and ContJJ(v,c) gate is a controlled gate with control qubit |u) and target 
qubit |c) that applies a single qubit unitary operator U on |c) only if |u) = |1), 
so every qubit |cf) carries a value of the constraint c; based on the values of Va 
and Vb in the balanced assignment \xk), i.e. the values of and respectively. 
Let \zk) = jcgCi ... c^_i), then the system can be re-written as follows. 
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(24) 


M-1 

l^'e) = a X] ® 

fc =0 

where every \xk) is entangled with the corresponding \zk)- The aim of the next 
stage is to find \zk) with the maximum number of |l)’s for the max-bisection 
problem or to find \zk) with the minimum number of |l)’s for the min-bisection 
problem. 


| 0 ) 

| 0 ) 

fo) 

| 0 ) 



Figure 3: Quantum circuits for (a)the MAX operator and (b) the MIN operator, 
followed by a partial measurement then a negation to reset the auxiliary qubit 
\aX2). 


3.3 Mciximization of the Satisfied Constraints 

Let I'l/'c) be a superposition on M states as follows, 

M-l 

\'^c) = Oi^\zk), ( 25 ) 

k=0 

where each \zk) is an m-qubit state and let dk = (zk) be the number of I’s in 
state \zk) such that \zk) ^ lO)®*”, i.e. dk ^ 0. This will be referred to as the 
1-distance of \zk)- 

The max-bisection graph \Xmax) is equivalent to find the state \zmax) with 
dmax = max{dk, 0 < k < M — 1} and the state \zmin) with dmin = min{dk, 0 < 
fc < M — 1} is equivalent to the min-bisection graph \xmin)- Finding the state 
\zmin) with the minimum number of I’s is equivalent to finding the state with 
the maximum number of O’s, so, to clear ambiguity, let dmaxi = dmax be the 
maximum number of I’s and dmaxo = dmin be the maximum number of O’s, 
where the number of O’s in \zk) will be referred to as the 0-distance of \zk)- 
To find either \zmax) or jzmin), when \ipc) is measured, add an auxiliary 
qubit \ax 2 ) initialized to state |0) to the system \tpc) as follows, 

\fpm) = li’c) 0 |0) 

M-l 

E kfe) ® |o) ■ 

fe =0 
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= a 


(26) 
































To main idea to find \zmax) is to apply partial negation on the state of \ax 2 ) 
entangled with \zk) based on the number of I’s in \zk), i.e. more I’s in \zk), 
gives more negation to the state of \ax 2 ) entangled with \zk)- If the number of 
I’s in \zk) is TO, then the entangled state of 10 x 2 ) will be fully negated. The 
TO*^ partial negation operator is the to‘^ root of ax and can be calculated using 
diagonalization as follows, 


V = 


1 

2 


1+t 1-t 

1-t 1+t 


(27) 


where t = ^—1, and applying V for d times on a qubit is equivalent to the 
operator. 


= 


1 

2 


1 + t’’* 1 - ■ 

1 +’ 


(28) 


such that if d = to, then = ax- To amplify the amplitude of the state 
\zmax), apply the operator MAX on lipm) as will be shown later, where MAX is 
an operator on to + 1 qubits register that applies V conditionally for to times 
on 10 x 2 ) based on the number of I’s in |coCi... Cm-i) as follows (as shown in 
Figure EKa)), 


MAX = Cont-V{co,ax2)Cont_V{ci,ax2) ■ ■ ■Cont_V{cm-i,ciX2), (29) 
so, if di is the number of c; = 1 in |coCi... Cm-i) then. 


MAX (|coCi...Cm_i) 0 |0)) = \cQCi...Cm-i) 0 { ^^2 2 

Amplifying the amplitude of the state \zmin) with the minimum number of 
I’s is equivalent to amplifying the amplitude of the state with the maximum 
number of O’s. To find \zmin), apply the operator MIN on |^m) as will be 
shown later, where MIN is an operator on to + I qubits register that applies V 
conditionally for to times on \ax 2 ) based on the number of O’s in |coCi... Cm-i) 
as follows (as shown in Figure EJb)), 

MIN = Cont_V (cq, ax 2 )Cont-V (c7, 0 x 2 ) • ■ • ContM (cm-i, < 1 x 2 ), (31) 

where cf is a temporary negation of c; before and after the application of 
Cont_V[ci,ax 2 ) as shown in Figure El so, if do is the number of c; = 0 in 
|coCi... Cm-i) then. 


M/A^(|coCi...Cm_i) 0 |0)) = |coCi...Cm_i) 0 f—^— |0) H-^— |I) j . (32) 

For the sake of simplicity and to avoid duplication, the operator Q will 
denote either the operator MAX or the operator MIN, d will denote either di 
or do, l^s) will denote either \zmax) or \zmin), and dg will denote either dmaxi 

or djnaxQ-, SO, 
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(33) 



-d 



Q (|coCi...Cm-i) 0 |0)) = |coCi...Cm-i) 0 


| 0 ) + 


and the probabilities of finding the auxiliary qubit \ax 2 ) in state |0) or |1) when 
measured is respectively as follows, 


Pr{ax2 = 0)= = cos^ (^) , 


(34) 


Pr{ax2 = 1) = ^ = sin^ {j£) . 


To find the state in \'\pm), the proposed algorithm is as follows, as shown 
in Figure [31 

1- Let \ij;r) = It/’m). 

2- Repeat the following steps for r times, 

i- Apply the operator Q on |'0r)- 

ii- Measure |ax 2 ), if |aa:: 2 ) = |1), then let the system post-measurement 
is \ipr), apply crx on \ax 2 ) to reset to |0) for the next iteration and 
then go to Step (i), otherwise restart the stage and go to Step (1). 

3- Measure the first m qubits in \'ipr) to read |2:s). 

For simplicity and without loss of generality, assume that a single jz^) exists 
in IV’d), although such states will exist in couples since each |zs) is entangled 
with a variable assignment jxs) and each jccs) is equivalent to |^), moreover, dif¬ 
ferent variable assignments might give rise to constraint vectors with maximum 
distance, but such information is not known in advance. 

Assuming that the algorithm finds \ax 2 ) = |1) for r times in a row, then the 
probability of finding |aa: 2 ) = |1) after Step (2-i) in the I'** iteration, i.e. r = 1 
is given by, 



(35) 


The probability of finding ji/jr) = l-^s) after Step (2-i) in the 1®* iteration, i.e. 
r = 1 is given by, 



(36) 


' \2m)' ^ ' 

The probability of finding |aa; 2 ) = |1) after Step (2-i) in the iteration, 
i.e. r > 1 is given by. 


M-l 


Pr^^\ax2 = 1 ) 



(37) 
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Figure 4: The probability of success for a max-bisection instance of the graph 
shown in Figure [I] with n = 8 and m = 12 where the probability of success 
of \ax 2 ) is 0.6091 after the first iteration and with probability of success of 
0.7939 after iterating the algorithm where the probability of success of \zmax) 
is amplified to reach the probability of success of \ax 2 )- 


The probability of finding \'ipr) = l-^s) after Step (2-i) in the iteration, 
i.e. r > 1 is given by, 

= z.) = — --, (38) 

E (fe) 

k=0 

To get the highest probability of success for Pr{'ijjr = Zg), Step (2) should 
be repeated until, \PA'^\ax 2 = 1) — PA^'^{tl}r = Zs)| < e for small e > 0 as 
shown in Figure SI This happens when ^/C5^«sin2’’ < e. Since the 

Sine function is a decreasing function then for sufficient large r, 


M-l 

E 

k—0,k^s 


Sin 


/ dfc7r \ 

V2my 


sin 


2r 



(39) 


where dns is the next maximum distance less than dg ■ The values of dg and dug 
are unknown in advance, so let dg = m he the number of edges, then in the 
worst case when dg = m, dns = m—l and m = n{n— l)/2, the required number 
of iterations r for e = 10“^ and A > 0 can be calculated using the formula. 


0 < sin^*' 


/ (m- l)7r \ 
\ 2m ) 


< e, 


(40) 
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( 41 ) 


r > 


_ log(e) 

2 log (sin 
_ log(lO-^) 
21og(cos(2^)) 

>A(2^)^ 

= O , 


)) 


where 0 < m < . For a complete graph where m = ; then the 

upper bound for the required number of iterations r is O (n^). Assuming that 
a single |zs) exists in the superposition will increase the required number of 
iterations, so it is important to notice here that the probability of success will 
not be over-cooked by increasing the required number of iteration r similar to 
the common amplitude amplification techniques. 


3.4 Adjustments on the Proposed Algorithm 

During the above discussion, two problems will arise during the implementation 
of the proposed algorithm. The first one is to finding \ax 2 ) = |1) for r times 
in a row which a critical issue in the success of the proposed algorithm to 
terminate in polynomial time. The second problem is that the value of dg is not 
known in advance, where the value of Pr^^\aX 2 = 1) shown in Eqn. 1351 olavs 
an important role in the success of Ending \aX 2 ) = |1) in the next iterations, 
this value depends heavily on the density of I’s, i.e. the ratio 

Consider the case of a complete graph with even number of vertices, where 
the number of egdes m = and all | 2 :fc)’s are equivalent and each can be 

taken as |zs) then, 


Pr^^\ax2 


1 ) 



(42) 


This case is an easy case where setting m = dg in rrd^ root of ax will lead 
to a probability of success of certainty after a single iteration. Assuming a blind 
approach where dg is not known, then this case represents the worst ratio ^ 
where the probability of success will be ~ 0.5 for sufficient large graph. Iterating 
the algorithm will not lead to any increase in the probability of both |zs) and 
10x2). 

In the following, adjustments on the proposed algorithm for the max-bisection 
and the min-bisection graph will be presented to overcome these problems, i.e. 
to be able to End |aa; 2 ) = |1) after the first iteration with the highest probability 
of success without a priori knowledge of dg. 


Adjustment for the Max-Bisection Problem 

In an arbitrary graph, the density of I’s will be jn the case of a complete 

graph, there are M states with 1-distance {dk) equals to This case represents 

2 

the worst density of I’s where the density will be 2 n{n-i) slightly greater than 
0.5 for arbitrary large n. Iterating the proposed algorithm will not amplify the 
amplitudes after arbitrary number of iterations. To overcome this problem, add 
Mmaa; temporary qubits initialized to state |1) to the register |coCi...Cr„_i) as 
follows, 
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Figure 5: The probability of success for a max-bisection instance of the graph 
shown in Figure [T] with n = 8, m = 12, fimax = 31 and <5 = 0.9, where the 
probability of success of \ax 2 ) is 0.9305 after the first iteration and with proba¬ 
bility of success of 0.9662 after iterating the algorithm where the probability of 
success of \zmax) Is amplified to reach the probability of success of \ax 2 )- 


ICqCi ...Cm_i) y IcqCi . . . CjYi—\C’fnCYYi-\.\ . . . , (43) 

so that the extended number of edges niext will be niext = m + Umax and 
V = will be used instead ot V = '^ox hr the MAX operator, then the 

density of I’s will be 2 n(n-i)+ 4 ii —• ^ probability of success Pxmax to 

find 1 0 x 2 ) = |1) after the first iteration, 

( '^ \ ”4 f Mmax ) \ 

PA^\ax2 = 1) = A/a^sin^ , ) ^ > Pw, (44) 

V2(^+fc«); 

then the required number of temporary qubits Umax is calculated as follows, 


Mmax > Y (2w - 1) - , (45) 

where uj = — sin 

TT 

P?'max = SMa'^ such that 0 < (5 < 1. For example, if i5 = 0.9, then Pr^^'> (ax 2 = 1) 
will be at least 90% as shown in Figure [51 To conclude, the problem of low den¬ 
sity of I’s can be solved with a polynomial increase in the number of qubits 
to get the solution \zmax) in O (wg^.^) = O (n^) iterations with arbitrary high 
probability J < 1 to terminate in poly-time, i.e. to read [ 0 x 2 ) = |1) for r times 
in a row. 




Mq2 


and Pxinax < with Ma^ = 1 so let 
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Adjustment for the Min-Bisection Problem 

Similar to the above approach, in an arbitrary graph, the density of O’s will be 
. In the case of a complete graph, there are M states with 0-distance (dfc) 

equals to ~ X- represents the worst density of O’s where the 

density will be slightly less than 0.5 for arbitrary large n. Iterating the 

proposed algorithm will not lead to any amplification after arbitrary number of 
iterations. To overcome this problem, add fimin temporary qubits initialized to 
state |0) to the register |coCi...Cm-i) as follows, 

IcqCi... 

Cm —1 ) -)■ |coCi . . . 

Cm —lCmCm+1 ■ • ■ — ^ ), (46) 

so that the extended number of edges rriext will be rriext = m + fimin and 
V = will be used instead of P = in the MIN operator, then the 

density of O’s will be 2 ra(w-i)+ 4 /r‘" • S^t a probability of success Pvmax to 
find \ax 2 ) = |1) after the first iteration, 


[ax2 


1) = Ma^ sin^ 



> Pr 


_ ' max? 


(47) 


then the required number of temporary qubits ^min is calculated as follows, 


n 

/^min ^ 


2a;- 1 
1 — a; 


n 

2 ’ 


(48) 


where u = f sin ^ and Pr^ax < Ma^, with Ma^ = 1 so let 

P?'max = 5Ma^ such that 0 < (5 < 1. For example, if i5 = 0.9, then Pr^^^ ( 0 x 2 = 1) 
will be at least 90%. To conclude similar to the case of the max-bisection graph, 
the problem of low density of O’s can be solved with a polynomial increase in the 
number of qubits, larger than the case of the max-bisection graph, to get the 
solution \zmin) in O {jn^xi^ — ^ iterations with arbitrary high probability 
(5 < 1 to terminate in poly-time, i.e. to read |ax 2 ) = |1) for r times in a row. 


4 Conclusion 

Given an undirected graph G with even number of vertices n and m unweighted 
edges. The paper proposed a BQP algorithm to solve the max-bisection problem 
and the min-bisection problem, where a general graph is considered for both 
problems. 

The proposed algorithm uses a representation of the two problems as a 
Boolean constraint satisfaction problem, where the set of edges of a graph 
are represented as a set of constraints. The algorithm is divided into three 
stages, the first stage prepares a superposition of all possible equally sized graph 
partitions in O ( -s/n) using an amplitude amplihcation technique that runs in 

O for A = 2" and M is the number of possible graph partitions. The 

algorithm, in the second stage, evaluates the set of constraints for all possible 
graph partitions. In the third stage, the algorithm amplifies the amplitudes of 
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the best graph bisection that achieves maximum/minimum satisfaction to the 
set of constraints using an amplitude amplification technique that applies an 
iterative partial negation where more negation is given to the set of constrains 
with more satisfied constrains and a partial measurement to amplify the set of 
constraints with more negation. The third stage runs in 0{m?) and in the worst 
case runs in O(n^) for a dense graph. It is shown that the proposed algorithm 
achieves an arbitrary high probability of success of 1 — e for small e > 0 using a 
polynomial increase in the space resources by adding dummy constraints with 
predefined values to give more negation to the best graph bisection. 
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